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Among the expected sources of gravitational waves for the Laser Interferometer Space Antenna 
(LISA) is the capture of solar-mass compact stars by massive black holes residing in galactic centers. 
We construct a simple model for such a capture, in which the compact star moves freely on a circular 
orbit in the equatorial plane of the massive black hole. We consider the gravitational waves emitted 
during the late stages of orbital evolution, shortly before the orbiting mass reaches the innermost 
stable circular orbit. We construct a simple model for the gravitational-wave signal, in which the 
phasing of the waves plays the dominant role. The signal's behavior depends on a number of 
parameters, including fi, the mass of the orbiting star, M, the mass of the central black hole, and 
J, the black hole's angular momentum. We calculate, using our simplified model, and in the limit 
of large signal-to-noise ratio, the accuracy with which these quantities can be estimated during 
a gravitational-wave measurement. For concreteness we consider a typical system consisting of a 
10 Mq black hole orbiting a nonrotating black hole of mass 10 6 Mq, whose gravitational waves 
are monitored during an entire year before the orbiting mass reaches the innermost stable circular 
orbit. Defining \ = cJ/GM 2 and rj = (i/M, we find Ax~5x l(T 2 /p, An/rj ~ 6 x lCT 2 /p, and 
AM/M ~ 2 x 10 _3 /p. Here, p denotes the signal-to-noise ratio associated with the signal and its 
measurement. That these uncertainties are all much smaller than 1/p, the signal-to-noise ratio level, 
is due to the large number of wave cycles received by the detector in the course of one year. These 
are the main results of this paper. Our simplified model also suggests a method for experimentally 
testing the strong-field predictions of general relativity. 
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I. INTRODUCTION AND SUMMARY 

A. Space-based gravitational-wave detectors 

The Laser Interferometer Space Antenna (LISA) 
project pj] was proposed in December 1993 as a corner- 
stone project for the "Horizon 2000 Plus" program of the 
European Space Agency. If the project is realized, LISA 
might start searching for gravitational waves in the year 
2010. 

The LISA interferometer consists of three pairs of 
spacecrafts, each pair located at the corners of an equilat- 
eral triangle. The constellation moves on a heliocentric 
orbit at 1 AU, following the Earth by an angle of 20°. 
Each side of the triangle has a length of 5 x 10 6 km, and 
the plane of the triangle is oriented at 60° with respect 
to the ecliptic. 

LISA uses laser interferometry to detect gravitational 
waves in the frequency band between 0.1 and 100 mHz. 
Such waves have much too low a frequency to be measur- 
able by Earth-based interferometers, which are severely 
limited by seismic noise at these frequencies LISA's 
projected sensitivity is such that a gravitational wave of 
amplitude 10~ 23 would give rise to a signal-to-noise ratio 
of 5 after one year of observation. 



Among the expected sources of gravitational waves rel- 
evant to LISA are solar oscillations || , binary systems of 
normal and compact stars within the galaxy Q, massive 
black-hole binaries in other galaxies B, and primordial 
gravitational waves Q . 

In this paper we consider another type of source, also 
expected to be highly relevant: the capture of solar-mass 
compact stars by massive black holes ||f7| . 

B. Stellar capture by massive black holes 

The most likely explanation for active galactic nuclei 
involves the tidal disruption of normal stars by a massive 
black hole residing in the galactic center. It is plausi- 
ble that compact stars (such as neutron stars and black 
holes) will also be captured by the massive black hole. 
However, such stars will not be disrupted, and their com- 
plex orbital motion will generate gravitational waves || . 
If the central black hole has a mass in the interval be- 
tween 10 5 and 10 7 solar masses, then these waves will 
have their frequency within the LISA bandwidth [[l] , and 
will be detected if they are sufficiently strong. 

It is expected that typically, a captured star will move 
on a highly eccentric orbit around the massive black hole 
|^-^| . Such orbits are easily perturbed by the other stars 
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in the galactic nucleus, and this can have a damaging 
effect on the gravitational-wave signal. The probability 
for an orbit to be perturbed decreases with increasing 
values of the orbiting mass. For this reason it is judged 
Jj|,[io| that the capture of compact stars with masses in 
the interval between 5 and 10 solar masses would be the 
most relevant for observations by LISA. Compact stars of 
such masses cannot be neutron stars, and would therefore 
be black holes. It was estimated [p]JTo|] that the rate of 
occurrence for such events could be as large as several 
per year, and that the signal-to-noise ratio could be near 
40 for a capture occurring at a distance corresponding to 
a cosmological redshift of one. 

In this paper we consider the gravitational waves pro- 
duced during the capture of a 10 Mq black hole by a mas- 
sive, 10 6 Mq black hole. We are specifically interested in 
what can be learned from observing these waves. 

C. Information carried by the waves 

The gravitational waves produced during the capture 
of a compact star by a massive black hole are rich in 
information. In principle, measuring the waves should 
allow us to extract the value of the black-hole parameters 
(mass and angular momentum), as well as to test the 
strong-field predictions of general relativity [0. In this 
paper we wish to propose a strategy for doing both. 

A possible method for extracting the black-hole param- 
eters, which we shall not explore in this paper, would go 
as follows g§. 

The gravitational waves generated by the bound mo- 
tion of a mass around a center come with a certain dis- 
tribution in frequency. For example, circular motion in 
the equatorial plane of a black hole produces waves with 
frequencies at every harmonic of the orbital frequency 
p3| . (The strongest waves come at twice the orbital fre- 
quency.) As an other example, the gravitational waves 
produced by a mass in eccentric motion (also in the equa- 
torial plane) come in harmonics of two fundamental fre- 
quencies @|||| , respectively associated with the angular 
and radial components of the motion. 

Precisely how the waves are distributed in frequency 
depends on the orbital parameters (for example, the ec- 
centricity). Measuring the distribution at some moment 
in time therefore leads to an estimation of these parame- 
ters. Now, because the system is losing energy and angu- 
lar momentum to gravitational waves, the orbital param- 
eters are not constant, but evolve in time. Precisely how 
they do so depends on the black-hole parameters (and 
also on the mass of the orbiting object). Measuring how 
the orbital parameters evolve in time therefore leads to 
an estimation of the black-hole parameters. 

Such a method was recently examined by Finn, Ori, 
and Thorne |T^ |, in the restricted context of circular, 
equatorial orbits. The method involves measuring the 
relative amplitude of each frequency component of the 



gravitational waves. In Sec. V we will see that amplitude 
measurements typically come with relative statistical un- 
certainties of order 1/p, where p is the signal-to-noise ra- 
tio associated with the signal and its measurement. (The 
uncertainty is due to the detector noise.) It follows that 
by using this method, the estimation of the black-hole 
parameters can be effected with a relative accuracy of 
order 1/p. 

Another method for extracting the black-hole parame- 
ters, also not explored in this paper, would involve mon- 
itoring that part of the gravitational-wave signal which 
corresponds to the black hole settling down after the or- 
biting mass' final plunge jl^Jl^]. This part of the sig- 
nal is strongly dominated by the black-hole quasi-normal 
modes fl6[] , whose frequencies and damping times are 
functions of the hole's mass and angular momentum. 
Such a method was examined by Echeverria 1T3] and Finn 
p"5| , who also conclude that the black- hole parameters 
can be estimated with a relative accuracy of order 1/p. 

The method considered in this paper is essentially due 
to Finn and Chernoff jTj]], and to Cutler and Flanagan 
p8[ . It goes as follows. 

As was pointed out above, the gravitational waves 
come in many frequency components. One of these, 
preferably the dominant one, is selected, and its phase 
is carefully monitored. Let us denote by / the frequency 
of this component. Because the system loses energy and 
angular momentum to gravitational waves, / increases 
slowly with time. The phasing of the waves is determined 
by df /dt, the rate of change of the frequency. In turn, 
this depends on the black-hole parameters (and also on 
the mass of the orbiting object). Monitoring the phase 
therefore leads to a way of estimating the black-hole pa- 
rameters. 

We will show in Sec. V that the estimation of those 
parameters which influence the phasing of the waves can 
be effected with a relative accuracy of order 1/ p c g , where 
the effective signal-to-noise ratio is given by p c g — 2ixM p. 
Here, N is the total number of wave cycles received by 
the gravitational-wave detector. This number will be es- 
timated in Sec. Ill: For a binary system consisting of 
two black holes, one of mass 10 Mq, the other of mass 
10 6 Mq, and for a signal monitored during an entire year, 
M ~ 10 5 . The method considered in this paper therefore 
leads to a much enhanced accuracy for the estimation of 
the black-hole parameters. 

D. Realistic captures and our simplified model 

The realistic description of a capture is complicated. 
The captured star moves on a highly eccentric orbit J^-||] 
which may not lie in the equatorial plane of the massive 
black hole. The orbital parameters (such as periastron 
and eccentricity) evolve in time, in part stochastically 
(due to perturbation by other stars) and in part secularly 
(due to gravitational radiation reaction). 
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To model a realistic capture, and to calculate the grav- 
itational waves emitted, is a difficult problem. To a first 
approximation, in which one ignores the external pertur- 
bation and the radiation reaction, the orbit is accurately 
represented by a bound geodesic of the Kerr spacetime 
||l9| , ^of . Such an orbit is characterized by three param- 
eters: orbital energy, orbital angular momentum, and a 
third called the Carter constant. 

An improved approximation would incorporate a small 
stochastic variation of these parameters, which would ac- 
count for the external perturbation, and also a small 
secular variation, which would account for the radiation 
reaction This approximation should be quite ac- 

curate if the density of stars is sufficiently low near the 
central black hole. 

The main difficulty preventing this problem from find- 
ing a solution resides in the fact that no prescription 
is currently available for calculating the secular evolu- 
tion of the Carter constant. (See, however, Refs. [^l] p3[ . 
The evolution of the other orbital parameters is obtained 
by computing how much energy and angular momen- 
tum are carried off by the gravitational waves. This is 
done by applying the methods of black-hole perturbation 
theory However, because of the system's small 

mass ratio, this problem should be much more tractable 
than the nonrestricted general relativistic two-body prob- 
lem, which currently is only amenable to post-Newtonian 
methods |J]. 

In this paper we choose to simplify the problem as 
much as possible, in order to introduce a simple, ana- 
lytic expression for the gravitational-wave signal. We 
shall assume that the captured object moves freely on 
a circular orbit in the equatorial plane of the massive 
black hole. We shall ignore all external perturbations, 
and focus on the secular evolution induced by the emis- 
sion of gravitational waves. Because the orbit is circular 
and equatorial, this evolution is dictated entirely by the 
loss of orbital energy. And because the ratio of the two 
masses is small, this loss is accurately computed using 
the methods of black-hole perturbation theory E3] . 



E. Model gravitational-wave signal 

As the system loses energy and angular momentum 
to gravitational waves, /, the frequency of the selected 
gravitational-wave component, slowly increases. For a 
circular orbit, the dominant component is at twice the 
orbital frequency; if f2 denotes the angular velocity of 
the orbiting mass, then 2irf = 20. It can be shown ||] 
that during most of the orbital evolution, the adiabatic 
approximation holds. This states that the radiation re- 
action timescale is much longer than the orbital period, 
so that df/dt <C l// 2 . Our analysis relies heavily on the 
adiabatic approximation. 

There exists an upper bound for the gravitational-wave 
frequency, which we denote by fo- This bound is set by 



the innermost stable circular orbit (ISCO) of the Kerr 
spacetime, and is equal to Qo/ir, where f2o is the an- 
gular velocity of a test mass moving on the ISCO. The 
gravitational waves therefore sweep up in frequency until 
/ = fo- At this point, the orbiting mass can no longer 
maintain a circular orbit, and is forced to plunge into the 
massive black hole, thereby shutting off the gravitational 
waves. It should be noted that the adiabatic approxima- 
tion necessarily breaks down as the orbiting mass reaches 
the ISCO @. 

We shall focus our attention on that part of the orbital 
evolution for which the gravitational waves are strongest 
and most interesting. This corresponds to the very late 
stages, shortly before the orbiting mass reaches the ISCO. 
If we assume that / is always very close to fo, then we 
may express all relevant quantities as Taylor expansions 
about the ISCO. 

We are mostly interested in the phasing of the waves, 
which is determined by df/dt. We will find in Sec. II B 
that this can be expressed as 



'1 
dt 



3 1 - f/fo 



l-2a(l-/// ) + 



(1.1) 



Here, V and a are dimensionless expansion coefficients 
which depend on /z, the mass of the orbiting star, M, 
the mass of the central black hole, and |x| = J/M 2 , its 
dimensionless angular-momentum parameter. (We use 
units such that G = c = 1; J denotes the black hole's spin 
angular momentum, and x will be defined more precisely 
in Sec. IV.) The dots denote terms of h igher order in 
1 — f/fo- It is to be noted that Eq. (LI) is singular at 
/ = fo- This reflects the breakdown of the adiabatic 
approximation at the ISCO. This issue is discussed in 
detail in Sec. II C. 

We use Eq. (.LI) to construct a model for the 



gravitational-wave signal originating from the capture of 
a compact star by a massive black hole. The model signal 
then depends on the parameters fo, V, and a, which cod- 
ify information about the source. We ask how well these 
parameters can be estimated during a gravitational-wave 
measurement, and what they can tell us about fi, M, and 
X- 

In the following we will refer to {fo, V, a} as the signal 
parameters, and to {/i,M,x} as the source parameters. 



F. A remark on our expansion approach 



The expansion (LI) is accurate only if 1 — f/fo is al- 
ways much smaller than unity. We must ask whether 
this is true for a system with = 10 M Q , M = 10 6 M , 
whose gravitational waves are monitored for an entire 
year before the orbiting mass reaches the innermost sta- 
ble circular orbit. The answer, most unfortunately, is 
no. For such a system, as we shall see in Sec. Ill, 
e = A///o ~ 1.1, where A/ is the frequency interval 
over which the signal is monitored. 
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To bu ild a model signal upon such expansions as 
Eq. ( LI ) therefore seems a bad idea if one is at all inter- 
ested in accuracy. Should we then give up this approach 
and do something better? We choose to give a negative 
answer: We will not give up, and we shall proceed with 
our expansions, even though the resulting signal may not 
be very accurate. 

The reason is that to provide a better representation 
for df/dt would require a great deal more work, and could 
not be done in the simple, analytic way advocated in this 
paper. Indeed, the gravitational-wave signal would have 
to be generated numerically. Now, numerical signals will 
have to be provided when an attempt is made to con- 
struct more sophisticated models, for example, models 
that incorporate highly eccentric orbits. However, we 
feel it would not be worthwhile to introduce numerical 
signals at this stage. 

Instead, for the sake of simplicity, we shall continue to 
treat 1 — ///o as a formally small quantity, and to carry 
out expansions in powers of that quantity. We hope that 
our conclusions will not depend strongly on this assump- 
tion, and that we can safely put e ~ 1 at the end of the 
calculation. 



G. The main results 

Our model gravitational-wave signal is characterized 
by the parameters fo, V, and a, which are functions of 
the source parameters p, M, and x- 

The signal parameters can be estimated during a 
gravitational- wave measurement. The statistical uncer- 
tainties associated with such an estimation are calculated 
in detail in Sec. V, in the limit of large signal-to-noise 
ratio. For concreteness we consider a typical system con- 
sisting of a 10 Mq black hole orbiting a nonrotating black 
hole of mass 10 6 Mq, whose gravitational waves are mon- 
itored for an entire year before the orbiting mass reaches 
the innermost stable circular orbit. For such a system we 
find that the measurement uncertainties are given by 



^ ~ 5.5 x 10- 4 /p, 
Jo 

— ~3.8x 10- 3 /p, 
Aa ~ 1.0 x 10~ 2 /p, 



(1.2) 



where p is the signal-to-noise ratio. A more precise, 
and more general, statement of these results appears in 
Eq. ( |5.12j ). These uncertainties should be compared with 
that associated with A, the signal's amplitude parame- 
ter: AA/A = 1/p. That the phase parameters can be 
estimated much more accurately is due to the large num- 
ber of wave cycles received by the detector. We indeed 
recall that TV ~ 10 5 . 



Ax ^ 5.0 x lQ- 2 /p, 
— ~ 5.5 x lQ-' 2 /p, 



(1.3) 



A more precise, and more general, statement of there 
results can be found in Sec. V C. We recall from Sec. I 
B that gravitational waves produced during a capture 
occurring at cosmological distances could be measured 
with p ~ 40. 



H. Testing general relativity 



Equations (1.2) imply that the source parameters \, 
rj = p/M , and M can be determined with uncertainties 



In order to obtain Eqs. fll.q ) from ( |l.2| ), the relation 
between the source parameters {p, M, \} and the signal 
parameters {fo,V,a} must be known. This relation is 
based upon two conditions: that the central mass is a 
Kerr black hole, and that the system emits gravitational 
waves according to the quantitative predictions of Ein- 
stein's theory. Indeed, general relativity is used to calcu- 
late dE/dt, the energy lost per unit time to gravitational 
waves. It is also used to write / = SI /it, the statement 
that the strongest waves come at twice the orbital fre- 
quency. On the other hand, the Kerr metric is used to 
calculate dE/dfl, the relation between orbital energy and 
angular velocity. Putting all this together gives df/dt in 
terms of the source parameters. 

It is easy to imagine a situation in which either one, or 
both, of these conditions would be violated. For exam- 
ple, the central mass could be something quite different 
from a Kerr black hole, or Einstein's theory could make 
wrong predictions regarding the emission of gravitational 
waves. As a specific example, we may recall that in the 
Brans-Dicke theory j2(| , the existence of scalar radiation 
implies a different result for dE/dt from what is pre- 
dicted by general relativity [^7f . However, by virtue of 
the uniqueness of the Kerr black hole pq ), the relation 
dE/dQ, is preserved in the Brans-Dicke theory. 

A breakdown of what we shall call the standard model 
(general relativity is valid and the central mass is a Kerr 
black hole) does not logically invalidate Eq. (1.1). In- 
deed, this equation is based solely on the assumptions 
that the orbiting mass moves on a circular orbit, and 
that the spacetime possesses an innermost stable circu- 
lar orbit. These conditions could easily be preserved in 
an alternative model. 

The standard model implies a particular relationship 
between the source and signal parameters. Another 
model would presumably imply a different relationship. 
We may therefore ask whether a measurement of the sig- 
nal parameters might allow us to test the validity of the 
standard model. This question is considered in Sec. V D, 
where we show that a measurement of a outside the in- 
terval (— oo, 1.8) would unambiguously signal the break- 
down of the standard model. 
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I. Conclusion 



II. MODEL GRAVITATIONAL- WAVE SIGNAL 



The analysis presented in this paper is by no means 
definitive, and needs to be improved in many respects. 
First, the crude representation (1.1) for df/dt will have 
to be replaced by something more reliable, presumably of 
numerical form. Second, the unrealistic assumption that 
the orbit must be equatorial and circular will have to be 
removed. All of this will require a great deal of effort. 

However, our simplified analysis should not be crit- 
icized too harshly. It is, after all, a first approach to 
a difficult problem. And our con clusions show promise: 
We believe that our results ( |l.3| ) for the measurement 
accuracies — most especially the fact that the phasing 
method yields uncertainties which are much smaller than 
1/ p — provide ample motivation for pursuing this work 
and adding new layers of complexity. 



J. Organization of this paper 

We begin in Sec. II with a detailed introduction of our 
model for the gravitational-wave signal. The assump- 
tions are reviewed and explained in Sec. II A. The signal 
is explicitly constructed in Sec. II B. Finally, some of the 
approximations employed are motivated in Sec. II C. 

In Sec. Ill we consider a specific case of the general 
framework introduced in Sec. II, in which the central 
mass is a Schwarzschild black hole. This section also 
provides numerical estimates for such quantities as fo, 
V, a, e, and J\f. 

In Sec. IV we generalize the discussion of Sec. Ill to 
the case of a Kerr black hole. 

Section V contains a discussion of the statistical un- 
certainties associated with the estimation of the signal 
parameters. In Sec. V A we re-introduce the model sig- 
nal. In Sec. V B, expressions are derived for the mea- 
surement uncertainties. These are converted, in Sec. V 
C, into uncertainties for the source parameters, assuming 
the validity of the standard model. Finally, in Sec. V D, 
a method is proposed to test the validity of the standard 
model. 

The Appendix contains a general discussion of orbital 
motion in spacetimes which are stationary, asymptoti- 
cally flat, axially symmetric, and reflection symmetric 
about the equatorial plane. (The orbits are restricted to 
this plane.) The equations of motion are introduced in 
Sec. A 1. Circular orbits are defined precisely in Sec. A 
2. In Sec. A 3 we consider two neighboring circular or- 
bits, and derive relationships between the displacements 
in radius, energy, angular momentum, and angular veloc- 
ity. Finally, the innermost stable circular orbit is defined 
precisely in Sec. A 4. 



In this section we introduce our model for the gravita- 
tional waves emitted during the late stages of orbital evo- 
lution of a binary system with small mass ratio, shortly 
before the smaller mass reaches the innermost stable cir- 
cular orbit associated with the larger mass. In Sec. II A 
we list and explain the assumptions involved in formu- 
lating the model. In Sec. II B we construct the signal. In 
Sec. II C we discuss the validity of some of the approxi- 
mations used in subsection B. 



A. Assumptions 

To formulate the model we shall make a number of 
assumptions. These are listed and explained below. 

Assumption 1: The binary system consists of a small 
mass fj, and a large mass M, so that fi/M <C 1. In typical 
situations we shall take /i = 10 M© and M = 10 6 M Q . 

Assumption 2: The spacetime geometry associated 
with the mass M is stationary, asymptotically flat, axi- 
ally symmetric, and reflection symmetric about the equa- 
torial plane. However, it is not assumed that M is a Kerr 
black hole, or even that the metric is a solution to the 
Einstein field equations. Although most of this paper will 
be devoted to the specific case of a Kerr black hole, we 
wish to keep the framework as general as possible in or- 
der to later identify a strategy for experimentally testing 
the predictions of general relativity. 

Assumption 2 leads to a number of simplifications. 
First, the fact that the spacetime is stationary and axi- 
ally symmetric implies that f2 = dc/)/dt, the angular ve- 
locity of the orbiting mass /it, is uniform if the orbit is 
circular. Here, 4> is the angle that goes around the axis 
of symmetry, and t is proper time for observers at rest 
at infinity. Second, the fact that the spacetime is reflec- 
tion symmetric about the equatorial plane implies that 
if started in this plane, the orbit will remain in the plane 
forever (whether it is circular or not). Consequently, we 
can further impose: 

Assumption 3: The mass fj, moves on a circular orbit 
in the equatorial plane of the mass M. This assumption 
is imposed solely for the sake of simplicity, and cannot 
be expected to hold in realistic situations. This was dis- 
cussed in detail in Sec. I D. 

Assumption 4-' The spacetime possesses an innermost 
stable circular orbit (ISCO). The ISCO is defined to be 
that particular circular orbit for which an infinitesimal 
change in orbital energy induces a large change in the 
orbital radius, or in the angular velocity Q. 

The ISCO can be thought of as the limiting member of 
a sequence of circular orbits, each labeled by the value of 
its angular velocity. For Q smaller than the critical value 
f^o, stable circular orbits exist; for f2 larger than Oo, 
circular orbits are no longer mechanically stable. Bound 
orbits within the ISCO all plunge toward the mass M. 
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Assumption 5: The binary system emits gravitational 
waves. It is not assumed that Einstein's theory gives 
the correct description of gravitational radiation, for the 
same reason as was expressed in the context of assump- 
tion 2. 

Due to the emission of gravitational waves, the system 
loses energy and angular momentum, and as a result, the 
mass fi undergoes orbital evolution: its angular velocity 
does not stay constant but keeps on increasing as the 
orbital radius decreases. 

Combining assumptions 1 and 5 leads to the conse- 
quence that while Q increases due to the radiation reac- 
tion, this increase occurs adiabatically ||: the timescale 
for the increase is much larger than the orbital period. 
The orbital evolution can therefore be represented by a 
succession of circular orbits (ordered by the increasing 
value of the angular velocity), with the mass fi moving 
quasi-statically from one orbit to the next. 

It is important to state that the adiabatic approxima- 
tion must break down as the orbiting mass reaches the 
innermost stable circular orbit Q . This follows from the 
very definition of the ISCO. Nevertheless, the adiabatic 
approximation holds for most of the orbital evolution, 
and its breakdown will not be an obstacle in the forth- 
coming analysis. 

Assumption 6: The gravitational waves emitted by the 
binary system are monitored by a space-based interfer- 
ometric detector (such as LISA) during an entire year 
before the mass /i reaches the innermost stable circular 
orbit. 

We shall calculate that during this time, the angular 
velocity changes by a fractional amount of order unity. 
We shall also find that Af, the total number of wave cycles 
received by the detector during the year's worth of ob- 
servation, is quite large. In typical situations, JV ~ 10 5 . 



B. Construction of the signal 

In brief, the model gravitational-wave signal is based 
upon a Taylor expansion about the innermost stable cir- 
cular orbit. The justification for this was given in Sec. I 
F, and will be repeated in Sec. II C. The model signal 
is also based upon the adiabatic approximation, which is 
itself a consequence of assumptions 1 and 5. And finally, 
we shall construct the signal by focusing mostly on its 
phasing, as was motivated in Sec. I C. 

We now proceed. Combining the adiabatic approxi- 
mation with assumption 3 (regarding the circularity of 
the orbit) implies that the gravitational waves have a 
frequency component at every harmonic of the angular 
velocity fl(t). (Here we indicate that the angular velocity 
evolves slowly with time.) More precisely pJlSj], the fre- 
quency F m associated with the m'th harmonic is given by 
2it F m (t) = mtt(t). We single out the strongest frequency 
component, and assume that only this component is ac- 
tually measured by the gravitational-wave detector. (In 



general relativity the dominant harmonic is at to = 2.) 
This particular frequency will be denoted by F(t). 

The phasing of the wave is determined by the function 
F(t). To obtain this we invoke once again the adiabatic 
approximation, and write 



dF 



dE/dt 

dt ~ dE/dF' 



(2.1) 



where E denotes the orbital energy of the mass fi. The 
numerator of Eq. (2.1) represents the energy lost per unit 
time to gravitational waves. The denominator gives the 
relationship between orbital energy and angular veloc- 
ity. The adiabatic approximation ensures that E and 
F change very little over the course of one orbit, so that 
these quantities are actually well defined. As was pointed 
out before, the adiabatic approximation, and therefore 
Eq. (2.1) also, break down as the mass \i reaches the in- 
nermost stable circular orbit. This does not represent an 
obstacle to our modeling; we will return to this point in 
Sec. II C. 

We denote by fo the gravitational-wave frequency at 
the innermost stable circular orbit. To indicate that fo 
must scale as the inverse power of M, we introduce the 
notation 



TtM/q 



(2.2) 



where vo is a scale-invariant quantity; the factor of n was 
inserted for convenience. 

We expand dE/dt about the innermost stable circular 
orbit. Since the gravitational-wave luminosity must scale 
as the square of the mass ratio, we write 



f = -(£)VMi-jki-f// o) + 



(2.3) 



Here, p and p are phenomenological parameters. The 
value that these parameters take depends on the nature 
of the mass M and on the correct description of grav- 
itational radiation. Specific cases will be considered in 
Sec. Ill and VI. For the time being, however, we leave 
these parameters undetermined. The factor of i>o 10 was 
inserted for convenience, and the dots designate terms of 
higher order in 1 — F/ fo. 

We also expand dE/dF about the innermost stable 
circular orbit. To do this we recall that the ISCO is 
defined to be that particular orbit for which an infinites- 
imal change in orbital energy produces a large change 
in the orbital radius, or in the angular velocity (which 
is related to the radius by a Kepler- type relation). The 
defining property of the ISCO is therefore that dE/dF 
must go to zero there. (See the Appendix for further de- 
tails, and for a justification of the fact that dE/dF must 
vanish linearly with 1 — F/ fo.) We also note that dE/dF 
must scale like the product of the two masses, as can be 
quickly seen from an elementary Newtonian calculation. 
We therefore write 

§ = q (1 - F/fo) [1 - q (1 - F/fo) + ■■■]. (2.4) 

db vo 



G 



As before, q and q are phenomenological parameters, and 
a factor of l/i>o was ins er ted for convenience. 
Combining Eqs. (0)-(EE|) yields 



^ = I V f 2 (1 - F/fo)- 1 [1 - 2a (1 - F/./o) + 



where 



3?rp 5 
q M '' 



V = — ^vo°, a = ±(p-q). 



(2.5) 



(2.6) 



Equation (|]^) can easily be integrated to give t(F), and 
the phase function 

$(F) = /" 2nF(t)dt = J 2nF-^dF (2.7) 

can be obtained similarly. 

We take the time-domain gravitational-wave signal to 

be 



h(F(t)) =A(l + --)e 



(2- 



Here, A is the signal's amplitude at the time the mass 
p, reaches the innermost stable circular orbit; A is in- 
versely proportional to the distance to the source. We 
shall also need an expression for the frequency-domain 
signal, which is obtained by Fourier transforming h(t): 

h(f) = J h{t) e 2 ™" dt = J h(F) e 2 ™/*^) J- dF. (2.9) 

The frequency / must be distinguished from F(t), the 
function of time which was also referred to as the fre- 
quency. 

We evaluate the Fourier transform by invoking the 
stationary-phase approximation p9| , which actually fol- 
lows from the adiabatic approximation |l8[ |. (We shal l 
return to this point shortly.) After substituting Eq. ( |2.g| ) 
into (2.S), we find that the argument of the exponen- 
tial becomes i times 2irft(F) — = f(F). Using 



Eq. (2/7), we see that this phase has a stationary point 
when = dip/dF = 2w(f - F)(dt/dF), that is, when 
F = f. The stationary-phase approximation then re- 
turns a factor [— v"(/)] -1 ^ 2 (1 — ///o) -1 ^ 2 which must 
be combined with the factor dt/dF\f oc (1 — f/fo) com- 
ing from inside the integral. After simple manipulations, 
including the absorption of constants into a redefinition 
of A, we arrive at 



h{f) =A(l + ■■■){!- f/fo) 1/2 e^\ 

V(/) = 27T/t - <l>0 + <t>(f), 
(1-///0) 3 



(2.10) 



V 



[1 + a(l ~ f/fo) + 



Here, the constants t = t(fo) and (po = $(/o) + 7r /4 
respectively represent time and phase at the innermost 



stable circular orbit. As before, the dots designate terms 
of higher order in 1 — /// p. 

We shall adopt Eq. ( |2.10| ) as our model for the 
(frequency-domain) gravitational- wave signal. It involves 
six phenomenological parameters: the amplitude A, the 
final time to, the final phase </>o, the final frequency fo, 
as well as V and a coming from the expansion of dF/dt 
about the innermost stable circular orbit. 



C. Validity of approximations 

Our model gravitational-wave signal is based upon a 
Taylor expansion about the innermost stable circular or- 
bit. Such an expansion can be accurate only if 1 — f/fo is 
small throughout the entire frequency interval over which 
the signal is monitored. If /, denotes the initial value of 
the frequency, we would require e = 1 — /j//o < 1. 

In Sec. Ill we will find that in typical situations, 
e ~ 1.1, and is therefore not small. This means that 
our model signal cannot be very accurate for frequencies 
near /j. This is most unfortunate, but we shall never- 
theless proceed with our model signal, which at least has 
the virtue of being extremely simple. In view of the fact 
that our model also incorporates unrealistic assumptions 
(that the orbit must be equatorial and circular), we feel 
that this degree of sophistication should be sufficient for 
our purposes. We do not expect that the use, in our cal- 
culations, of a more exact representation for dF/dt would 
lead to very different conclusions. A more complete dis- 
cussion of this point was given in Sec. I F. 

We now discuss the breakdown of the adiabatic and 
stationary-phase approximations. The adiabatic approx- 
imation becomes invalid when the timescale over which 
F(t) changes becomes comparable to the orbital period. 
This translates to (1/F) (dF/dt) ~ F or, upon using 
Eq. (ph and F ~ /„, 



(1 - F/fo) ~ V. 



(2.11) 



It can be shown that the same criterion applies to the 
breakdown of the stationary-phase approximation Jl8| ] . 

To circumvent the problems associated with the break- 
dow n of these approximations, we shall impose that 
Eq. fl2.10D is valid for 



f<f f = fo (1 - V) 



(2.12) 



only. For / > //, the signal is taken to be identically zero . 
In any event, the part of the signal for which Eq. ( 2.12 ) 
is violated turns out to be irrelevant for the purpose of 
Sec. V. Indeed, we shall see that most of the signal-to- 
noise ratio is accumulated when ( 1 — f/fo) is much larger 
than V. This is reflected in Eq. (^J), below. 
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III. SCHWARZSCHILD BLACK HOLE 

To give concreteness to the framework introduced in 
the preceding section, and to prepare the way for the 
next, here we consider the special case in which the mass 
M is a Schwarzschild black hole. In this and the following 
section we assume that general relativity gives the correct 
description of gravitational radiation. 

In this context, the signal's dominant frequency is 
at twice the orbital frequency, so that 2irF — 2f2 = 
2(M/r 3 ) 1 / 2 , where r is the orbital radius in Schwarzschild 
coordinates. Since the inner mos t stable circular orbit is 
located at r = 6M H, Eq. (O) gives 



- fi-1/2 



v = 6 



0.4083. 



(3.1) 



The gravitational waves emitted by a mass fj, in circu- 
lar motion near the ISCO of a Schwarzschild black hole 
of mass M can be calculated using the methods of black- 
hole perturbation theory (l^j3^,^l| . The numerical eval- 
uation of dE/dt then yi elds estimates for the constants p 
and p appearing in Eq. ( |2.3| ) . Using the results presented 
in Ref. |32(1, we obtain 



p~ 7.30, 



P 



3.78. 



(3.2) 



On the other hand, dE/dF can be calculated analyt- 
ically, since this quantity relates orbital energy and an- 
gular velocity for circular geodesies of the Schwarzschild 
spacetime. An elementary calculation yields 



E = n{l - 2u 2 )(l - 3v 



(3.3) 



where v = (M/r) 1 ' 2 = (Mil) 1 / 3 = (ttMF) 1 / 3 . After 
differentiation, 



dE 
dF 



TTflM 

3w 



(l -6v 2 )(l -3v 2 ) 



-3/2 



(3.4) 



Fina lly, expansion about F = fo = (6 3 / 2 7rM) 1 yields 
Eq. (UM, with 



q = ~ 1.9746, 



q= 2- 



(3.5) 



We now use the results obtained thus far to estimate 
the parameters which characterize the gravitational-wave 
signal. 

We b egin with the final-frequency parameter, fo- Use 
of Eqs. (^2) and (3A) yields the relationship between f 
and M, which we write as 



fo = 4.40 



t'o 



M 



mHz. 



(3.6) 



v 0.4083/ V 1Q6M 0. 

We have chosen a fiducial value of 1 x 10 6 M Q for M, to 
ensure that for such a mass, f = 4.4 x 10 -3 Hz, where 
LISA's sensit ivit y is near maximum. In the way it is 
written, Eq. (3.6) remains valid even when the mass M 



is not a Schwarzschild black hole; in such a situation vq 
would differ from 6~ 1//2 . The same will be true for most 
of the equations liste d b elow. 



Next, we use Eqs. (gj), (J3JJ), (^2|), and <^§) to obtain 



V = 3.95 x 10" 



M 



10 6 Mr. 



7.30 

-l 



''o 



1.9746 

5 



10 Mp 



0.4083 



and 



1.62. 



(3.7) 



(3.8) 



Continuing with our estimates, we consider a situation 
in which the gravitational waves are monitored during 
an entire year before the mass fj, reaches the innermost 
stable circular orbit. We ask: what is the correspond- 
ing fractional change in frequency? To answer this we 
integrate Eq. (O) for t(F): 



t(F) = t. 



2^V/ 



(l-F//o) 2 (1 + •••). 



(3.9) 



Choosing fi to be the initial frequency, letting At = to — 
t(fi), e = 1 — fi/ fo, and ignoring the corrections terms 
in Eq. (pi]), we obtain 



: = ( y V/oAi 



1/2 



= 1.07 







4.40 mHz 



1/2 



V 



3.95 x 10" 



1/2 



At 

yr 



(3.10) 

1/2 



We recall that according to Eqs. (2.12) and (3.7), the 
signal must be cut off at a final frequency // given by 



1 - ff/fo = V ~ 10- 



(3.11) 



Finally, we estimate Af, the number of wave cycles re- 
ceived by the detector during the year's worth of obser- 
vation. Even though the signal is not monochromatic, 
for our purposes it is sufficient to define 



M = f At 



1.39 x 10 5 



3e 2 



fo 



4.40 mHz 



At 

yr 



(3.12) 



2ttV 



The last line follows from Eq. ( 3.10 ) an d is di splayed here 
for future reference. We see from Eq. (3.12) that for an 
integration time of one year, the total number of wave 
cycles is quite appreciable. 
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IV. KERR BLACK HOLE 

We now generalize the discussion of the previous sec- 
tion to the case of a Kerr black hole. 

As before we have that the dominant component of the 
gravitational waves has a frequency given by F — 0,/tt. 
For a Kerr black hole, this reads [EOl 



nMF = 



1 



.3/2 



X 



(4.1) 



Here, we use a dimensionless radial coordinate r which 
is a rescaled version of the usual Boyer-Lindquist coor- 
dinate, namely, r = r^/M . The dimensionless constant 
X is related to J, the intrinsic (spin) angular momentum 
of the black hole, and also to L, the unit vector pointing 
in the direction of orbital angular momentum: 



-0.5 



0.5 




0.4 




0.3 




0.2 




0.1 





0.5 

X 



FIG. 1. A plot of 7rM/o = t>o 3 as a function of \- 



X 



J L 

M 2 ' 



(4.2) 



Thus, x is positive if the orbital motion proceeds in the 
same direction as the black-hole rotation (such an orbit 
will be termed direct); x is negative if the orbital mo- 
tion proceeds in the direction opposite to the black-hole 
rotation (such an orbit will be termed retrograde). It 
should be noted that x is restricted to the interval [—1,1]; 
|x| = l describes an extreme Kerr black hole. 

The innermost stable circular orbit is located at the 
value of r which satisfies r 2 — 6r + %x r1 ^ 2 ~ 3x 2 = 0- The 
solution is [E0| 



3 + B 



X 



A=l+(1- X 



V(3- A)(3 + A + 2B), 

(i + x) 1/3 + (i-x) 1/3 



2\l/3 



(4.3) 



B= v%C 2 



A 2 . 



We observe that for x — 0, f = 6, the Schwarzschild 
value. For x = — 1 (extreme black hole, retrograde orbit) 
we have r — 9, while f — 1 for x = 1 (extreme black 
hole, direct orbits). The equation r = 1 also describes 
the event horizon of an extreme Kerr black hole. It is 
useful to record the expansion 



f=l + 2 2 / 3 (l- X ) 1/3 + ^(l-x) 2/3 

+ Ig (i-x) + -, 

which is valid for x ap p roachi ng unity. 
Combining Eqs. ([2.2]) and (4.1) gives 



(4.4) 



7rM/ EE Vq 3 



1 



F3/2 



■ X 



(4.5) 



A plot of vq 3 (x) is provided in Fig. 1. We note that for 
X = 0, vq 3 reduces to the Schwarzschild value, 6 -3 / 2 ~ 
0.0680. For X = ~1, vq 3 = 1/26 ~ 0.0385, while v 3 = 



1/2 for x = 1- We also calculate, for orbits approaching 
the innermost stable circular orbit, 



1 - F/fo 



r 3/2 _ f 3/2 
r 3/2 _|_ 
3f l/2 



2(f 3 / 2 



(r 



3(5f 3 / 2 - X ) 

(r — r) 



8r 1 / 2 (f 3 / 2 +x) 
This series can be inverted to give 

2(r 3 / 2 + X ) , 



(4.6) 



3f!/2 

(5f 3 / 2 



(1 - F/fo) 
x) (f 3/2 + ; 



9f 2 



(1-Wof 



(4.7) 



We will use this result in order to reproduce, for the 
speci fic c ase of a Kerr black hole, the expansions ( |2.3| ) 
and (2.4) for dE/dt and dE/dF, respectively. 

The computation of the gravitational- wave luminosity, 
_Egw = —dE/dt, for a specified circular orbit of the Kerr 
spacetime, is carried out numerically using the methods 
of black-hole perturbation theory ^5|. In such calcula- 
tions, the value of the radius is conveniently used to char- 
acterize the orbit. By computing £gw(^) and i?Q W (f), 
where a prime denotes differentiation with respect to r, 
we have access to the parameters p and p. The precise 
relations arc 



P= — 



EGw(r) 



vo 1 



(4.8) 



E' GW (f) 



3fi/2 t)Q 3 EGW (r) 







TABLE I. The tabulated form of p and p as functions of 
X- These values were obtained from numerical data provided 
by Masaru Shibata. 



x 


P 


V 


-0.9 


7.79 


3.76 


-0.8 


7.75 


3.77 


-0.7 


7.70 


3.78 


-0.6 


7.65 


3.78 


-0.5 


7.60 


3.78 


-0.4 


7.55 


3.78 


-0.3 


7.49 


3.78 


-0.2 


7.43 


3.78 


-0.1 


7.36 


3.78 


0.0 


7.30 


3.78 


0.1 


7.21 


3.77 


0.2 


7.12 


3.80 


0.3 


7.02 


3.77 


0.4 


6.90 


3.77 


0.5 


6.75 


3.76 


0.6 


6.57 


3.75 


0.7 


6.33 


3.72 


0.8 


5.95 


3.66 


0.9 


5.22 


3.46 



These parameters are listed in Table I for several values 
of X- The data from which p and p were obtained was 
kindly provided by Masaru Shibata (Osaka University, 
Japan). Because it is increasingly difficult to calculate 
i?Gw for increasing values of |%|, such computations were 
not carried out for \x\ > 0.9. 

As for the case of a Schwarzschild black hole, the cal- 
culation of dE/dF, and its expansion near the innermost 
stable circular orbit, proceeds analytically. 

We begin with the relation between orbital energy and 
radius l|20fl: 



r 3/2 _ 2r l/2 + x 
r 3/4( r 3/2 _ 3r l/2 + 2x) 1/2 ' 



(4.9) 



After differentiation and upon using Eq. (4.1), we obtain 

6r + 8xr l/2_ 3x 2 )(r 3/2 +x) 2 



dE _ nfiM (r 2 
dF ~ 3 



r 9/4( r 3/2 _ 3r l/2 + 2x) 3/2 



(4.10) 



We recognize, in the numerator, the expression r 2 — 6r + 
g^ r i/2 _ 2^2 w j 1 j cn vanishes when r — f. An exp ansion 
about this point will therefore be of the form (f2~4), as we 
shall see presently. 

We point out that this statement is false in the case of 
direct orbits around an extreme Kerr black hole. When 
X = 1 the denominator of dE/dF also goes to zero at 
r = f = 1, and as a consequence, dE/dF does not vanish 
at r = f. Instead, 




x=i, 



9\/3 



(4.11) 




FIG. 2. Plots of q (higher curve) and q (lower curve) as 
functions of X - Both quantities diverge as x — * 1- 

We may therefore conclude that the condition r = f does 
not, in the case of direct orbits around an extreme Kerr 
black hole, correspond to a genuine innermost stably cir- 
cular orbit, since dE/dF fails to go to zero there. The 
case x = 1 therefore represents a singular limit in the 
framework introduced in Sec. II, and we shall exclude 
it from further considerations. More details about this 
point can be found in the Appendix. 



As was stated before, expa nsio n of Eq. (4.10) about 
r — f and substitution of Eq. (4/7) reproduces Eq. (^4|), 
with 



4?r 



g Uo 8 f 13/4 ( f 3/2 _ 3f l/2 + 2x ) 



1/2 ' 



and 



15f 2 + ixr 1 ~ 39xr i/2 + 22x 

6f 3/2( f 3/2_ 3f l/2 + 2x ) 



(4.12) 



(4.13) 



Plots of q(x) and q(x) are provid ed in Fig. 2. We ob- 
serve that for x = 0, Eqs. ( 4.1 2| ) and ( |4 .13 ) red uce to 
the Schwarzschild expressions given by Eq. (3.5). For 
X = -1, q ^ 1.6401 and q = 17/81 ~ 0.2100. For 
X approaching ■ 
q ~ (2 7 /3/ 3 ) (1 



-1, q - (16-\/l27r/27) (1 - x)~ 1/3 and 
view of Eq. ([4.1l|), it is 



not surprising that these quantities should diverge in this 
limit. 

We use our results for vq, p , p, q, and q to calculate V 
and a according to Eq. ( |2.6[ ). Plots of these quantities 
are provided in Figs. 3 and 4. 

Be fore le aving t his se ction we point out that Eqs. 
( |3.7| ), and (3.1C)-( 3.12| ), which provide estimates for /q, 
V> /i//oj ff/fo, and Af, are quite general: they therefore 
apply equally well to the case of a rotating black hole. 



V. ESTIMATION OF SIGNAL PARAMETERS 
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FIG. 3. A plot of (M /fx) V as a function of %■ The com- 
puted data points at \ = {~ 0.9, —0.8, ■ • • ,0.9} are joined by 
straight line segments. 
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FIG. 4. A plot of a as a function of X- The computed data 
points at x = {—0.9, —0.8, • • ■ , 0.9} are joined by straight line 
segments. 



A. Introduction 

The last two sections were devoted to special cases of 
the general framework introduced in Sec. II. In such cases 
there exists a specific relationship between the signal pa- 
rameters {/o, V, a} and the source parameters {p, M, x}- 
For example, in the case of a Kerr black hole, a is a func- 
tion of the single variable \; this function is represented 
in Fig. 4. 

We now return to the general framework of Sec. II, in 
which the nature of the mass M is left unspecified. The 
framework also makes no specific choice regarding the 
correct description of gravitational radiation. 

In this framework the gravitational- wave signal is given 
by Eq. ( |2.10D , which we rewrite as 



h(f) = A(l-f/f ) 1/2 e^f\ 



^(/) = 2irft Q - 0o + *(/), 
(W//o) 3 



(5.1) 



V 



l + a(l-/// 



For the purpose of this section we consider Eq. (|5.1| ) to be 
exact, ignoring all relative corrections of order (1 — ///o) 
in the amplitude, and all relative corrections of order 
(1 — ///o) 2 in the phase. 

The signal is characterized by six parameters, which 
we group into the vector 



6 = (In .A, In / , In V, a, t , <f> ). 



(5.2) 



The purpose of this section is to calculate, for a given 
signal strength, how accurately the parameters 6 can be 
estimated during a gravitational- wave measurement. 



B. Parameter estimation 

The statistical theory of signal detection and measure- 
ment fUJH] nas repeatedly been applied to the specific 
case of gravitational waves |^,^6). For a complete ex- 
position, we refer the reader to Refs. [ p^|JT^j37[| . That 
the theory must be statistical in character follows from 
the facts that gravitational- wave detectors are noisy, and 
that the noise is a random process, here assumed to be 
stationary and Gaussian. In this subsection we apply the 
theory to calculate, in the limit of large signal-to-noise 
ratio, the expected statistical errors associated with the 
estimation of the signal parameters 6. 

We first introduce the inner product (-|-): If a and b 
are two functions of time, with Fourier transforms a and 
b, then 



(a\b) 



~a*(f)b(f)+a(f)b*(f) 

Sn(f) 



df. 



(5.3) 



Here, an asterisk denotes complex conjugation, and 
S n (f) is the spectral density of the detector noise, the 
Fourier transform of the noise's autocorrelation function. 
The detector has good sensitivity where S n (f) is small, 
and poor sensitivity where S n (f) is large. For LISA, 
S n (f) is small in the frequency bandwidth between 1CP 4 
and 10" 1 Hz §. The precise behavior, and overall nor- 
malization, of S n (f) will be of no concern to us. 

The signal-to-noise ratio p associated with the mea- 
surement of a signal h(t) by a detector with spectral 
density S n (f) is given by p 2 — (h\h). This is simply eval- 
uated for our model signal (^j]) . Expanding the spectral 
density as S n {f) = S n (f )[l+O(l- f / / )], neglecting the 
second term for consistency, and integrating from / = /j 
to / = //, we obtain 



2/o |-4| 



S n (fo 



(5.4) 



Here, e = 1 — fi/fo, and is given b y Eq (3.10). We have 
neglected V 2 in front of e 2 , cf. Eq. ( |3.11 ). From Eq. (5.4) 
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we see that the signal-to-noise ratio is proportional to the 
signal's amplitude times the square root of its duration 
(as measured by e), and inversely proportional to the 
noise's amplitude. 

We calculate the anticipated accuracy with which the 
parameters 8 can be estimated during a gravitational- 
wave measurement. In the limit of large signal-to-noise 
ratio (or small statistical errors), the calculation involves 
computing, and then inverting, the Fisher information 
matrix T a i, — (h ta \h t b), where h <a — dh/dO a is the partial 
derivative of the signal with respect to t he p arameter 9 a . 

With the ordering provided by Eq. ( ]5.2| ), the partial 
derivatives are given by 



where there is no summation over the repeated indices. 

To invert the Fisher matrix we proc eed a s follows. Af- 
ter computing the matrix we use Eq. (3.12) to express V 
in terms of e and Af. The matrix is then inverted alge- 
braically, using the symbolic manipulator Maple. Thus 
far all expressions are exact, and quite large. To sim- 
plify, we take advantage of the fact that according to 
Eq. ( 3.12 ), Af is numerically very large. We therefore 
expand our exact expression for S in powers of l/Af, 
keeping only the leading-order term. It can be verified 
that the resulting matrix is equal to Sigma, apart from 
a fractional correction of order l/Af 2 . 

We obtain the following results: 



where 



.9i = 1, 



92 

93 
.94 



1 - X 
2x 



g a h, 



2ix 3 , 
1 + — (3 + Aax) 



(5.5) 



V 



(1 + ax), 



IX 



.95 = i(27r/ )(l - x), 
.96 = -i, 



(5.6) 



and x = 1 — f/fo- Then, using Eq. (5.4) and the same 
approximation for S n (f) as was used before, we have that 

2 r-e 

r a h = — / [g* a (x)g b {x)+g a {x)gl{x)]xdx. (5.7) 
£ Jv 

These integrations are easily carried out. Because its 
expression is rather large, the Fisher matrix will not be 
displayed here. 

The meaning of the Fisher matrix is as follows |I^|L8| ■ 
Because of random noise, the output of a gravitational- 
wave detector is also random, even when a signal of the 
form h(t; 6) is known to be present. Consequently, an es- 
timation of the parameters 6 comes with inevitable sta- 
tistical uncertainty. It can be shown that in the limit 
of large signal-to-noise ratio, the probability distribu- 
tion function for the parameters 6 is of the Gaussian 
form exp[-ir ah (6» a - 8 a )(6 b - 6 b )}, where denotes the 
true value of the signal parameters. It follows that the 
variance-covariance matrix, E ab = ((6 a - 6 a )(8 b - 6 b )}, 
where (■) denotes an averaging with respect to the spec- 
ified distribution, is given simply by the inverse of the 
Fisher matrix, 



(r 



-l\ab 



(5.8) 



In terms of this, the "one-sigma" statistical error associ- 
ated with the estimation of parameter 9 a is given by 



A9 a = 



(5.9) 



AA 

A 
A/o 

fo 
AV 

V 

Aa 



1 

P 

280\/3 
2nAfp ' 

1890^ 
2irAfep ' 

378^ 



At = 



2vrAAe 2 p' 
105^2 
2irf ep ' 
105V2c 



(5.10) 



where 



a(e, a) = 1 



32 128 
21 ~ "63~' 
16 128 
27 ~ "8T Q 



£ 

256 
243 



a 2 e\ 



b(e, a) = l + yae+^(2 + 17a)a e 2 

+ — (1 + Aa)a 2 e 3 + — (1 + 4a) V e\ 



4 1 9 

c(e) = 1 eH e . 

w 21 98 



(5.11) 



To put this in a more concrete form we e valuat e a, b , 
and c at e = 1.1 and a = 1.6; cf. Eqs. (3J3) and ( 3.10 ). 
This gives a ~ 3.8, b ~ 170, and c ~ 0.80. We then 
choose Af = 1.4 x 10 5 and fo = 4.4 mHz as fiducial 
values. Substituting all this, we obtain 



AA _ 1 
~A~~~P 

Af _ 5.5 x 10" 4 
fo P 
AV _ 3.8 x 10" 3 
"V" ~ ~P 



Af 



1.4 x 10 5 

1/2 



3. 



e 

n 



Af 



Aa 



1.0 x 10" 2 



1/2 



170/ V 1-1 



1.4 x 10 5 

Af ' 
1.4 x 10 5 



12 



At 
A0 O 



81 min 
P 

19 cycles 
P 



a VYjl 

4.4 mHz / V 1 - 1 

\ 1/2 / \ -1 
_C_\ 

0.80/ V 1.1 



(5.12) 



From these expressions we observe that: (i) The rel- 
ative accuracy with which A, the amplitude parameter, 
can be estimated is precisely equal to the inverse signal- 
to-noise ratio. We shall call this level of accuracy the 
"signal-to-noise ratio level", (ii) The estimation of fo 
benefits the most from the large number of wave cycles; 
its accuracy is typically three orders of magnitude better 
than the signal-to-noise ratio level, (iii) The estimation 
of V and a also benefits from the large value of AT; their 
accuracies are respectively three and two orders of mag- 
nitude better than the signal-to-noise ratio level, (iv) 
The uncertainties associated with the parameters to and 
</>o are independent of M. 



C. Measurement of black- hole parameters 

For the purpose of this subsection we return to the 
special case in which the mass M is a Kerr black hole 
(Sec. IV). The question we ask is: What does a measure- 
ment of the signal parameters {fo, V, a} tell us about the 
source parameters {fi,M, x}? 

Inspection of Figs. 1, 3, and 4 reveals that the quan- 
tities nMfo = vo 3 , ?/ -1 V, and a are all functions of the 
single variable x- Here we have defined r\ to be the mass 
ratio: 



V = T7 



fj_ 

w 



We therefore write 



ttM/o = A( x ), 
V - 1 V = B( X ), 
a = C( X ). 



(5.13) 



(5.14) 



It follows from these relations that a measurement of a 
yields the value of x- Knowing x> a measurement of V 
then gives rj. And finally, a measurement of fo returns 
the value of M. In this way, the source parameters can 
all be estimated from a measurement of the gravitational- 
wave signal. The question remains as to how accurately 
all this can be done. This is the question we turn to next. 

We first consider the estimation of x, t he b lack-hole 
rotation parameter. It follows from Eq. (5.14) that in 
the limit of small measurement uncertainties, the error 
in x is given by 



1 



\C'(X)\ 



Aa, 



(5.15) 



where a prime denotes differentiation with respect to \. 
To express this in a more concrete form it would be de- 
sirable to explicitly compute the function C"(x)- Unfor- 
tunately, this cannot be done reliably, since C(x) was 
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FIG. 5. A plot of A' /tyA as a function of \- 



evaluated at only a few widely separated points. Never- 
theless, we crudely estimate C"(x) to be 0.2 at the some- 
what representative point X = 0. This yields 



Ax^5.0Aa forx = 0. 



(5.16) 



It can be seen from Fig. 4 that the uncertainty in x would 
be larger than this for negative values of x, while it would 
be smaller (and going to zer o for x — ► 1) f° r positive val- 
ues. Using Eq. fTT^ ), Eq. (E3|) means that for typical 
situations, Ax — 5.0 x 10~' z /p. The rotation parame- 
ter can therefore be estimated with a degree of accuracy 
approximately two orders of magnitude better than the 
signal-to-noise ratio level. 

We now repeat this procedure for rj, the mass-ratio 
parameter. In the limit of small uncertainties, Eq. (5.14) 
implies 



Arj 
V 



B'(X) 



B(X) 



Ax +ir - 



(5.17) 



A crude estimate reveals that \B'/B\ ~ 1.1 at x = 0. 
This means that the uncertainty in x dominates the un- 
certainty in t), and we find 



An 

— '- ~ l.lAx for x = 0. 
V 



(5.18) 



For typical situations, Ar\jr\ ~ 5.5 x 10~ 2 /p, a degree of 
accuracy one order of magnitude better than the signal- 
to-noise ratio level. 

Finally, we consider the accuracy with which M, the 
black-hole mass, can be estimated. We once more use 



Eq. (5.14) to derive 



AM 



Ax- 



A/q 
fo ' 



(5.19) 



A plot of A'/irA is provided in Fig. 5. Once more we 
find that the uncertainty in x dominates over the other 
contribution. Evaluation at x = yields 
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— ~ 0.238A X for X = 0. (5.20) 

It can be seen from Fig. 5 that the uncertainty in M 
would be smaller than this for negative values of \, 
larger for positive values, and significantly larger for \ 
approaching unity. For typical situations, AM/M ~ 
2.4 x 10~ 3 /p, a degree of accuracy three orders of mag- 
nitude better than the signal-to-noise ratio level. 

The bottom line is that during a gravitational-wave 
measurement, the source parameters {/x,Af, x} can be 
estimated with a degree of accuracy that is much better 
than the signal-to-noise ratio level. 

D. Testing general relativity 

Since the framework of Sec. II makes no explicit ref- 
erence to a particular model for the mass M, nor to a 
particular theory of gravitation, we may ask whether it 
leads to a viable strategy for testing the strong-field pre- 
dictions of general relativity. More precisely, we ask: Can 
a measurement of the signal parameters {fo, V, a} allow 
us to test the hypothesis according to which the central 
object is a Kerr black hole and the correct theory of grav- 
itation is general relativity? 

Clearly, measuring fo and V if of no use in this purpose. 
This is because fo and V respectively depend on M and 
r\ = IJ>/M, which are not known prior to the measure- 
ment. Therefore, measuring fo and V offers no way of 
distinguishing two different models for the source which, 
by adjusting the values of M and /i, give identical values 
for fo and V. 

The parameter a is more useful. It is a consequence 
of the standard model (the mass M is a Kerr black hole 
and the correct theory of gravitation is general relativ- 
ity) that a is restricted to the interval I sm = (— oo, 1.8). 
Measuring a therefore leads to a way of testing general 
relativity, since a measured value of a outside I sm would 
unambiguously invalidate the standard model, 

Of course, while viable in principle, this test is very 
much a partial one, and leaves much to be desired. For 
example, it is impossible to say which aspect of the stan- 
dard model is violated when a £ I sm . Another shortcom- 
ing is that there could well exist a wide class of alternative 
models for which a is restricted to an interval I which is a 
subset of I snl . The test offers no way of ruling out models 
belonging to this class. Finally, we must also recall that 
the entire framework is based upon a rather restrictive 
assumption: that the orbit must be circular and within 
the equatorial plane. 
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APPENDIX: 

In this Appendix we consider the motion of a test mass 
in a spacetime satisfying the following properties: it is 
stationary and asymptotically flat, possesses axial sym- 
metry, and is reflection symmetric about the equatorial 
plane. We do not assume that the spacetime metric satis- 
fies the Einstein field equations. However, for simplicity 
we shall assume that the test mass is moving freely, so 
that it follows a geodesic of the spacetime. Generaliza- 
tion to the presence of external forces (which also satisfy 
the symmetry requirements) would be straightforward. 
We shall focus our attention on equatorial, circular or- 
bits, the meaning of which will be made precise below. 

The results obtained in this Appendix form the basis 
for the construction of the model signal of Sec. II. 

1. Equations of motion 

We introduce coordinates (t, r, 8, <ft) to chart the space- 
time. We take t to be the time coordinate associated with 
the timelike Killing vector £ = d/dt, and normalized such 
that it is equal to proper time for observers at rest at 
infinity. We take r to be a well-behaved radial coordi- 
nate, such that circular orbits (defined in a coordinate- 
invariant way below) are described by r = constant. We 
take 8 and cj) to be the usual spherical coordinates, such 
that <f> is associated with the Killing vector x — d/d<fi. In 
these coordinates, according to the specified symmetries, 
the only off-diagonal component of the metric is g^. For 
the specific case of a Kerr black hole, the coordinates em- 
ployed reduce to the Boyer-Lindquist coordinates |HJ. 

Reflection symmetry about the equatorial plane en- 
sures that the motion of a test mass can be restricted to 
lie in this plane. (Of course, such a restriction represents 
a loss of generality). We therefore take 8 = ir/2, and the 
four-velocity is u a = dx a /dr — (t, f, 0, 0), where an over- 
dot denotes differentiation with respect to proper time r. 
The quantities 

E = -u£, L = u-x (Al) 

are constants of the motion, respectively representing or- 
bital energy per unit rest mass, and orbital angular mo- 
mentum per unit rest mass. 

It is easy to show that Eqs. ( |Al[ ) imply 
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t = T(r, E, L) = [g^E + g t4> L)/D, 
4> = <t>(r,E,L) = - {gt+E + g tt l) /D, 



(A2) 



where D = —gttg<j><j> + 9t4> 2 - On the other hand, the nor- 
malization condition u ■ u = — 1 gives 



= R(r, E,L) = 

Qrr 



1 g u E 2 + 2g t4> EL + g tt L 2 



g rr D 



(A3) 



As was indicated, the metric functions are considered 
to be functions of r only, since they are evaluated at 
= 7r/2. 

We define a function Q. by 



eft T(r, E, L) 



(A4) 



This represents the angular velocity of the test mass, as 
measured by an observer at rest at infinity. 



dR 
dE 



SE 



dR 

di 



SLn = 0. 



This eq uation takes a more familiar form if we use 
Eq. (A3) to evaluate the partial derivatives. We quickly 



arrive at 



5Eo = Qq 5Lq 



(A6) 



Equation ( |A5| ) also implies that d r R{rx,Ei,L\) = 0. 
Expanding, using Eq. ( |A5| ) to evaluate the partial deriva- 
tives, and Eqs. (A4) and (A6) to simplify the result, we 
arrive at 



d 2 R 



Or 2 



6r 



Qq g r 



d<P 8T 
"n 

or or 



SEn 



(A7) 



We notice that {d r <$>-Vld r T)/T = d r fl. Regularity of the 
radial coordinate therefore implies that the right-hand 
side of Eq. (A7) is nonzero whenever T ^ 0, or g rr =/= oo. 

A similar calculation would yield <5f2o m terms of either 
one of <5ro, SE , or SL . 



2. Circular orbits 



4. Innermost stable circular orbit 



We define an orbit to be circular if r is constant along 
the orbit. The condit ion s for this to be true are r = r = 
0. According to Eq. (A3), this means 



R(r,E,L) = — (r,E,L) = 0. 
Or 



(A5) 



These equations determine the values that E and L must 
take in order to produce a circular orbit at radius r. For 
an orbit at r — ro, we shall denote these values by Eq 
and Lq. 

Let £Iq be the angular-velocity function evaluated at 
r = ro, E = E (r ), and L = Lo(tq). We therefore 
have the important result that the angular velocity is uni- 
form along a circular orbit. This provides a coordinate- 
invariant definition of circularity, a property that has 
nothing to do with the geometric appearance of the or- 
bit. As a condition on our choice of radial coordinate, we 
require that Qq be a smooth, monotonic function of ro . 



3. Neighboring circular orbits 

We now consider two neighboring circular orbits, the 
first at r — ro (with associated parameters Eq and Lq), 
the second at r = ri (with associated parameters E\ and 
Li). We assume that n = ro + <5ro, with <5ro a small 
quantity. Correspondingly, E\ = Eq + SEq and L\ = 
Lq + SLq. We want to calculate various relationships 
between the displacements 5r , 5E , and 5L - 

From Eq. (A5) we must have R(r±, Ei, Lx) — 0. Ex- 
panding this about (ro, Eq, Lq), up to first order in the 
displacements, gives 



We assume that the spacetime is such that stable circu- 
lar orbits can only exist at radii larger than some limiting 



value 



At 



circular orbits become unstable, 



in the following sense: an infinitesimal change in the or- 
bital energy (or orbital angular momentum) produces a 
finite change in the orbital radius. The innermost stable 
circular orbit is therefore defined to be that particular 
circular orbit at which 



SEq 

Sr 

We shall assume that the condition 
U g rr \dr dr ) . 



(A8) 



(A9) 



holds. [For the specific case of a non-extreme Kerr black 
hole, this condition can be shown to hold for both direct 
and retrograde orbits, for any value of the hole's angular 
momentum. However, for the special case of an extreme 
Kerr black hole, the condition is violated for direct orbits. 
This is because <? rr |is C o = 00 j which signals the presence 
of an event horizon at r = ri sco . As a result, SE/dro 
never goes to zero for direct circular orbits around an 
extreme Kerr black hole.] 

Combining Eqs. (A8) and (|A9D, we obtain 



8 2 R 



dr 2 



= 0. 



(A10) 



The innermost stable circular orbit corresponds to an 
inflection point of the effective potential R. 
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We have required that f^o be a smooth, monotonic 
function of tq . Equation ( A8 ) therefore also implies that 



SEo 

sn 



= 0. 



(All) 



This is a coordinate-invariant characterization of the in- 
nermost stable circular orbit. 

It is easy to check that near ro = Hsco, 8Eq/8Q,q oc 
d 2 R/dr 2 \o. We assume that d 3 R/dr 3 does not vanish at 
To = fisco, so that the right-hand side vanishes linearly 
with ro — ri SC0 . We then quickly arrive at the following 
estimate: 



SEo 

sn 



k(flo — 



(A12) 



for some con stant k. This result provides the justification 
for Eq. (O), above. 
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